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Abstract. The structural inversion asymmetry- induced spin-orbit interaction of 
conduction band electrons in zinc-blende and wurtzite semiconductor structures is 
analysed allowing for a three-dimensional (3D) character of the external electric field 
and variation of the chemical composition. The interaction, taking into account all 
remote bands perturbatively, is presented with two contributions; a heterointerface 
term and a term caused by the external electric field. They have generally comparable 
strength and can be written in a unified manner only for 2D systems, where they 
can partially cancel each other. For quantum wires and dots composed of wurtzite 
semiconductors new terms appear, absent in zinc-blende structures, which acquire the 
standard Rashba form in 2D systems. 



PACS numbers; 73.21.-b, 73.40.Kp. 
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1. Introduction 

Systems with lowered space symmetry are generally characterized by spin-split energy 
states. This is a manifestation of the relativistic interaction of moving magnetic angular 
momentum with electric field [1]. In the physics of semiconductor nanostructures, we 
conventionally identify two reasons for the effect: lack of inversion symmetry of the unit 
cells of the constituent materials [21 [3l HJ [5] and the presence of the structural inversion 
asymmetry (SIA) on macroscopic scale, much larger than the unit cell [El [7]. 

Major results in this field have been obtained by using only symmetry arguments, 
with the method of invariants [8], which has become the most efficient instrument to 
studying electron states in semiconductors [9]. However, this method has a drawback, 
for its phenomenological nature does not yield a distinct link between bulk materials 
and heterostructures. Each reduction in dimensions, and more generally, each loss of 
elements of symmetry including a result of the application of strong enough external 
fields, creates a new system that calls for an independent analysis [IDl [HI 112] ■ 

Another successful phenomenological tool, the k ■ p method, the only extra 
requirement of which is mean- field approximation [9], does provide the missing link 
along with its limits of applicability [13]. For two-dimensional (2D) electron systems, 
this method has already proved to be capable of specifying the SIA mechanisms of the 
spin-orbit interaction and identifying them as the heterointerface induced and external 
electric field induced [TJ [m [151 [HI [H] • Knowledge of their strengths is crucial to band 
structure engineering aimed at manipulation with the spin degree of freedom, which 
finds important applications [HI [191 120] ■ 

An interesting situation holds for quantum dots. Their lateral confinement is 
usually modelled as parabolic, validating the inclusion of only ID interface spin-orbit 
interaction, see, e.g., [21]. However, attempts to make allowance also for the spin-orbit 
interaction due to a lateral field have been made in the past . This may be decisive 
for quantum dots based on piezoelectric materials such as GaN/AlN, where complex 
distributions of strain generate 3D pictures of strong (of the order of 10^ V/ cm) 'external' 
electric fields [231 1211 [23 (211 I2Z]- In addition to structures with 3D confinement, 
even nominally ID and 2D electron systems have their realizations in arrangements 
with boundaries and contacts inducing electric fields acting on the interior of the 
systems. That is why a ID spin-orbit interaction model may fail there. Indeed, it may 
underestimate by an order of magnitude [28] the electron spin splitting in GaN/AlN 
quantum wells, see e.g. |29j. 




We use the k ■ p method to systematically include the effect of an external electric 
field of arbitrary profile and general 3D variation in chemical composition on the spin- 
orbit interaction in zinc-blende and wurtzite heterostructures. Previous studies dealt 
with either 3D composition profile and arbitrary external electric field, considering them 
in detail only for cubic materials [H], or with ID electric fields for wurtzite materials 
[28] ■ In the former work, the basis functions comprised spin-orbit interaction. As a 
consequence, two different effects, describing the position-dependent effective mass and 
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governing the spin-orbit interaction, combined, and some of the spin-orbit terms were 
not singled out. An incomplete picture of the spin-orbit interaction resulted in the 
conclusion that the heterointerface contribution is small if the spin-orbit interaction 
parameter weakly varies with composition. We find it more convenient to use the 
original effective-mass method's spinless basis |30] to obtain all the necessary terms. 
After analysis of them, we found that the heterointerface contribution should always be 
taken into account. 

The study of wurtzite [28], borrowing the method used for zinc-blende materials 
[TJ [HI [17], was limited to the eight-band Kane-type model [31]. In that method, 
small (valence-band) envelope functions are excluded, resulting in a conduction 
band Hamiltonian that parametrically depends on its own eigenenergy [32]. Apart 
from neglecting remote bands whose contributions have never been evaluated 
(with questionable relevance for wide-bandgap materials), such a Hamiltonian has 
applicability problems, e.g. when time-dependent external fields are considered. Our 
resulting SIA spin-orbit interaction terms for conduction band states near the Brillouin 
zone centre in zinc-blende and wurtzite semiconductors are Hermitian and energy- 
independent. Any remote bands can be taken into account if proper material parameters 
are known. 

This paper is organised as follows. In section 2, we make a perspective analysis 
of the heterointerface- and external electric field-induced SIA mechanisms and find 
that they should generally be considered on equal footing. In section 3, we introduce 
a multi-band system of envelope-function equations used in section 4 to derive spin- 
orbit terms entering conduction band envelope-function equations, with details being 
given in appendix A. The general expressions for the spin-orbit interaction terms are 
then analysed for zinc-blende and wurtzite heterostructures. We discuss the results in 
section 5. 

2. Heterointerface- and external electric field- induced terms: a comparison 

The heterointerface spin-orbit contribution first appears in the third order jH] of the 
Lowdin perturbation scheme [9l|33]. It is proportional to the difference in the spin-orbit 
interaction parameters for the semiconductors of the heteroj unction. Analysing the 
'exact' expression for the parameter of the interface term in the Kane model [3 [T6l [T7] . 
we see that there should also be a contribution due to the valence band offset, present 
even when the spin-orbit interaction parameter does not vary with composition. Such a 
term, proportional to the band offset and the spin-orbit interaction parameter, will be 
available only if we consider the fourth order of the perturbation scheme. 

The electric field-induced contribution, which is proportional to the external electric 
field and the spin-orbit interaction parameter, arises only in the fourth order of the 
perturbation scheme if the Kohn-Luttinger basis functions pD] do not include spin-orbit 
interaction. Otherwise, Leibler showed that inclusion of the spin-orbit interaction in 
the basis makes both contributions present as third-order corrections [H]. For typical 
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Figure 1. A single-heterojunction potential well, showing the ground eigenenergy and 
its envelope function. 



heterostructures, the spin-orbit interaction energy is less or of the order of the band 
offsets. This, together with conventional arguments discussed in section 4, makes the 
original method's spinless basis [30] preferable. Then, such a perturbative classification 
of the interface and electric field-induced terms might indicate that the former term is 
stronger than the latter. 

To find out whether this is so, let us consider them for a quantum well grown 
along the 2;-axis and having for simplicity only one heteroj unction located at 2 = (see 
figure 1). With only Rashba-type spin-orbit contributions present, the effective mass 
equation for the electron envelope functions which correspond to the steady state 
eigenenergies e„, can be written as H2D^n = ^n^m with the Hamiltonian: 

H2D = ^ + U{z) + W{z) + (^R,6 (z) + ^2^^) {Kay - kya,) , (1) 

where is the momentum operator, m* is the effective mass, U{z) = Q{—z)6Us and 
W{z) are the potential energy of on electron in the crystalline potential and external 
scalar potential, respectively, Q{z) is the Heaviside step function, 6Us is the conduction 
band offset, 6{z) is the Dirac delta function, cr = (cr^, ay,az) is the Pauli matrix vector 
and Ri and R2 are material parameters. The term proportional to -Ri (i?i-term) defines 
the interface spin-orbit interaction. It is small by the parameter presented with a sum 
SA/Eg-\-A6Uy/ Eg, where ^ A is the difference of the valence band spin-orbit splittings for 
the two materials forming the heterojunction, A is the valence band spin-orbit splitting 
for the quantum well material, SUy is the valence band offset and Eg is the bandgap. The 
term proportional to R2 (-R2-term) is the electric field-induced spin-orbit interaction. It 
is small by a parameter that is of the order of A{W) / Eg, where (W) is a characteristic 
potential energy of an electron in the external scalar potential. 

Following [T7], we consider the Hamiltonian H2D0 of the zeroth-order approximation 
in the spin-orbit interaction: 

H2Do = ^ + Uiz) + Wiz). (2) 
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The corresponding envelope functions are designated as $on = |^): -f^2_Do^on = ^n^on- 
Using the eigenfunctions \n) of the Hamiltonian H2DQ and detaihng diagonal matrix 
elements of the commutator [A;2 , -^200] j which are zero, we immediately obtain 



holding for any subband n. Comparing this identity with the spin-orbit interaction 
part of ([1]), we see that the terms proportional to -Ri and R2 can both be written as 
a single interface contribution, found by Pfeffer and Zawadzki [17]. We argue that 
the identity (|3]) suggests that both spin-orbit terms can also be written as a single one 
proportional to the external electric field. This makes them mutually indistinguishable if 
phenomena involving intersubband transitions are not of interest, so that only a single 
effective material parameter can be extracted from experiment. In our opinion, the 
form requiring knowledge of the average electric field is preferable. The electric field 
can easily be estimated from the electrostatics of the semiconductor system involved, 
but it is impossible to make a direct evaluation of the magnitude of the envelope function 
at the heterointerface without numerical simulations. Analysis of a two-heterointerface 
quantum well does not alter the above conclusion. 

The identity ([3]) also helps us to make a simple comparison of the Ri- and i?2-terms, 
which involves the analysis of two small indicating parameters: 6 A/ Eg + ASU^/Eg and 
A6Us/Eg, respectively. As typically \SUs\ ~ \SUy\, we immediately conclude that if the 
spin-orbit interaction energies for the materials of the structure differ essentially, that 
is \SA\ ~ A, the heterointerface term dominates, as found by Leibler [H]. Otherwise, 
if |5A| ^ A, both terms should generally be retained, contrary to the conclusion in 
[H]. As examples, consider two popular semiconductor pairs forming heterostructures: 
GaAs/AlAs and GaN/AlN, taking the material parameters from [31]. For GaAs/AlAs, 
A ^ 0.34 eV, 5 A ^ -60 meV, Eg ^ 1.52 eV, 5Us ^ 1 eV and 5Uy ^ -0.53 eV. For the 
pair GaN/AlN, A ^ 14 meV, 6A ^ 5 meV, Eg ^ 3.5 eV, 6Us ^ 2 eV and 6U^ ^ -0.7 eV. 
The indicating parameters are comparable, |5A| ^ A, so we should keep both terms. 
Below, in section 4, we produce a more accurate evaluation for these two semiconductor 
pairs and show that both systems have an 'accidental' set of parameters leading to 
somewhat different conclusions. In the above estimates, we do not take into account 
remote bands whose effect has never been evaluated. For electron states in quantum 
wires and dots, the -Ri- and i?2-type terms cannot be reduced to a unified form because 
the operators and/or ky do not commute with the Hamiltonian, and the identity ([3]) 
cannot be directly used in the spin-orbit interaction terms of the Hamiltonian. 

In summary, the heterointerface- and external electric field-induced terms can be 
of comparable strength and, while being mutually indistinguishable for 2D electron 
systems, they cannot be written in a unified form for quantum wires and quantum dots. 
These make them be discrete and equally important. 




(3) 
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3. The multi-band system of envelope- function equations 

We now proceed to derive the SIA-governed spin-orbit terms using the k ■ p method. We 
consider a heterostructure composed of two semiconductors with the net single-particle 
potential energy of electron, which will be called 'potential' for brevity, U = U (r): 

U = Ui + f[U2- Ui] = Ui + f6U, (4) 

where Ui = Ui (r) and U2 = U2 (r) are the periodic lattice potentials of the nominally 
potential well material and the barrier material, respectively. We suppose that the 
phenomenological function / = /(r), which defines the profile of the structure, is of 
the order of unity or less even at the heterointerfaces [TU |35]. Ideally, / can even be a 
step-like function taking the values / = in the region of the potential well material, 
and / = 1 in the region of the barrier material [13]. It can have a variation in ID, 
2D or 3D to represent a potential well, a quantum wire or a quantum dot, respectively. 
The final expression will have a local character allowing plain generalization for an 
arbitrary composition. We set the only requirement that the semiconductors composing 
the structure are not too dissimilar, so that 5U can be treated as a small perturbation 
as compared to the basis potential Ui. 

In the mean-field approximation, the Schrodinger equation with the relativistic 
spin-orbit interaction term is [1]: 

+ U+ , + W]^{v) = e^> (r) . (5) 



2mo Arn^c^ 

Here mo is the free electron mass and c is the velocity of light in vacuum. The external 
scalar potential W = W {r) is weak, and we neglect its direct relativistic effect. 

Dealing with states near the Brillouin zone centre, it is convenient to use the 
complete set of Kohn-Luttinger functions [30j, with the unit-cell normalized periodic 
parts M„o = (r) = \n) specified as 
f fi^k^ \ 

+ Ui \n) = eno\n), (6) 



V 2mo 

where e„o is the nth band edge energy. The relativistic effect of the potential Ui is not 
included in the basis, but processed as a perturbation. We define the nth band envelope 
function in r-representation as An = An (r) with 

vl/(r) = ^AJn), (7) 

n 

where the summation is over all bands. 

Consider the following k ■ p system in r-representation, which is obtained by 
treating the functions / (r) and W {r) as 'gentle' [30]. This is a standard procedure 
of the envelope-function method [35], which neglects all 'central cell'-like corrections 
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due to a rapid variation in the function / (r) at the heterointerfaces |13j . 



Uno + 7^ + W r U„ + V ^An. 

V 2mo / ^ mo 

n' 



(8) 



n' 

where 




(9) 



(10) 



We define the matrix elements: p„„/ = (n | /ik | n'), 5f/„„/ = (n | 5f/ | n') and 
[Vf/i X k]^^, = {n I [Vf/i X k] | n'). The all-band system of equations ([8]) is valid for 
slowly varying envelope functions An [9]. 

We do not explicitly include the strain Hamiltonian into consideration, important 
for lattice- mismatched pairs, in particular for GaN/AlN. The proper procedure, detailed, 
e.g., in |36] (see also [37]), would lead to redundant complications, not essential for our 
results. It suffices to take into account here that the piezoelectric field due to strain 
along with the possible spontaneous polarization field has contributed to the 'external' 
potential W (r). For attainable values of strain, the infiuence of the deformation 
potentials on the band edge energies e^o, 'offsets' 5f/„„/ and the matrix elements -f^^lf"'* 
and n'f^f' is too weak to be included in the spin-orbit terms being derived. 

Ignoring the central cell corrections, some of them contributing to the spin-orbit 
interaction [13], we should be aware that there is a number of heterointerface-related 
effects that cannot be accounted for [12]. The corresponding material parameters are 
generally not expressed via bulk parameters of the constituent materials, and depend 
on microscopic structure of the heterointerface and its crystalline orientation [38] . 
Their estimates are scarce [39l HQl SI], and evidence that they produce a noticeable 
contribution to the spin-orbit interaction is currently absent. Possible speculations that 
they alone could explain the huge electron spin splitting in GaN/AlN quantum wells 
[29] have yet to have some grounds in first-principle band-structure calculations. 

We also omitted the /c- linear spin-orbit interaction term due to the potential t/i, 
see [21112]. In the third perturbation order (with two operators ^kp„„//mo), it generates 
/c-cubic bulk inversion asymmetry spin-orbit interaction terms, which are conventionally 
called the Dresselhaus term in zinc-blend materials and the Rashba term for wurtzite 
[43] . Their expressions are known [2[|11]. One could allow for the position dependence 
of these terms originating from a fc-linear spin-orbit interaction contribution due to 
the potential f{r)6U, which is also omitted from equation ([8]). Its inclusion may be 
consistent only for special cases of weakly localized electron states, where the barrier 
penetration is very essential. This term describes bulk inversion asymmetry position- 
dependent spin-orbit interaction. Also, we will not consider fc-linear bulk and the related 
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position-dependent terms for wurtzite structures appearing through the second-order 
perturbation (due to H^^"^ -\- f (r) H^^°^ and ^kp„„//mo) on the same grounds as above: 
the expression for the bulk term is known, see |42j, and an analogous position-dependent 
term has a very week effect on strongly localized states. 



4. Single-band spin-orbit Hamiltonian 

To allow for the interface spin-orbit interaction, we should deal with at least the third- 
order perturbation term that is a 'product' of f (r) H^^^"^ and two /ikp„„//mo. For 
a simple conduction band with the index m = s, the correction of the 3rd order in 



perturbation H' is, see Appendix A, equation (1A.14I) : 



^i? = ^{H2,S2h., (11) 

where the braces stand for a commutator, and H2 is a part of the perturbation H' having 
only non-diagonal couplings of the s-band with remote bands; S2 is given by flA.16p . 
Using i?if*°'* = Pss = 0, which holds both for zinc-blende and wurtzite materials at the 
r point of the Brillouin zone, with the help of expression (IA.16I) . we have: 

H(^) = \2EkEkEi±, (12) 

^ ^sl ^sV 

To obtain the SIA spin-orbit term originating due to the external potential 
W{v) supplemented with the crystalline potential f{r)6Unn', we deal with the fourth 
perturbation order in a product of H^^"\ W{v)6nn' + / (i") ^Unn' and two /ikp„„//mo. 
Here 5„„/ is the Kronecker delta. We do not consider a contribution from / (r) H^^f* as 
put in the place of H^^"\ It may produce a noticeable effect only if two conditions are 
satisfied: we treat weakly localized states with high enough probability for the electron 
to be in the barrier material, and H^^f^ ~ ^nn"^ ■ '^^^ latter condition would mean that 
the interface contribution (fT2|) dominated over the electric field-induced one (in terms of 
section 2, 6 A ~ A) canceling the very consideration of anything farther than the third 
perturbative order. 

If we included the spin-orbit interaction H*^^*°^ in the zero-order Hamiltonian Hq, 
we would not have to go down to the fourth perturbative order, as only the third will be 
required [H], expressed as a product of W{r) + / (r) 6Unn' and two ^kp„„//mo. Then, 
to be consistent, we would have to use the spinor basis functions that could not be 
treated as zero-order combinations of the functions Uno [30]. Instead, they should be 
constructed using not less than the first-order functions: 

<0 = «nO + V^^U„,o, (13) 



where the summation does not include degenerate terms. Using the functions f[T3|) with 
the third-order perturbative expressions [H], the final result will still have a character 
of the fourth-order smallness. We consistently treat H*^^*"-' as a perturbation to conform 
with the conventional classification of the basis functions. 
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The correction of the fourth order is, see appendix A, equation (lA.14p . 

Hif = \ {H2, S3},, - ^ {{{H2, Si} , Si} , Si},, , (14) 

where Si and S3 are given by ( 1A.15I) and (1A.17I) . Again, using the properties 
Hsl^°^ = Pss = and the expressions (lA.lSP and flA.17p . we have 



(15) 



^(4) ^ H'^iHwH[,i„H[„^ - i ( + ^ 

Another property we will use below is that the matrix element 5Unn' is finite only for 
functions \n) and \n') of the same symmetry. 

4-1. Zinc-blende semiconductors 

For the conduction band in structures based on zinc-blende semiconductors, expressions 
f[T^ and f[T^ are greatly simplified as we can neglect the matrix elements H^j"^"^ and 
jy^^so) prove it, we use the following arguments. Firstly, the spin-orbit 

interaction operator can be written as the product of the spin and orbital angular 
momentum near the atomic nuclei, where the interaction is essential. Secondly, the 
zone-centre function u,o, which transforms accordingly to the Fi representation of the 
space group T^, is composed of spherically symmetric atomic s orbitals with zero angular 
momentum. 

The functions u„o can be chosen real, so that 



SUnn' = SUn'n, Pnn' = " Pn'n, H^^, = — -f^„/„ . (16) 

Using these properties and changing the band summation indices, we have the following 
for the third perturbative order from ( fT2l) : 

^s>-i:%g^iv./*.-v./M. (17) 

where a, P = x, y, z, and summation over these indices is implied here and henceforth. 

For the fourth perturbative order, only the first term on the right-hand side of 
equation (fTS!) gives a finite contribution: 

ivl" ^"^0 

Finally, we arrive at the traditionally looking spin-orbit interaction terms [Ij, now 
generated by both the external scalar potential and variation in the chemical composition 
of the structure, in the form first given in pUE], Hgo = Hsf + H^f: 

Hso = RizB [V/ (r) X k] ■ + R2ZB [^W (r) x k] ■ o", (19) 
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where 

^^4 



rilZB — > A- A 2 

2im^ c^ujsi ujsv ujsi" 



(20) 



and 



p ^ fi\P.)sl ([V^l X k] j„, (Py),, ( 1 ^ 1 ^ 

Note that the hetero interface contribution, which is proportional to the parameter Rizb, 
originates not only due to the difference in the spin-orbit interaction energies, as given 
with the first term in the right-hand side of ( 120|) . but also due to the finite matrix 
elements SUui, see the second term of (1201) . All items are present in [H], but the latter, 
because of the chosen spinor basis, contributed to the position-dependent effective mass 
and was not analysed. 

For illustrative purposes, let us limit ourselves to the truncated eight-band Kane 
model, with the degenerate valence band edge's functions uxo, uyo and uzo transforming 
as X, y and z, respectively, in accordance with the Fis representation of the space group 
Trf. If we introduce P = i {s \ hkx \ X), Eg = Usx, SU^ = 6Uxx and 



3i Atuq (?■ ' 3i 4mQ (? 



(22) 



we obtain 



(23) 



and 



'imlEl V ' Ey 

R2ZB — 7. 2 173 ' (24) 

which coincide with the known result (see [71 [161 lH] and put the eigenenergy entering the 
Hamiltonians there e = e^o); for ^so, If^Asol, \5Uv\ -C Eg and ID external electric field. 
Using GaAs/AlAs band parameters [31], see also section 2, we have Rizb ~ 4 eVA^ 
and R2ZB ~ —6 A^. Note the different signs of these parameters. They partially cancel 
each other for GaAs/AlAs quantum wells, which is seen if we use the identity Q. If 
the conduction band offset were 5Us = —Rizb/R2ZB ~ 0.7 eV, they would cancel each 
other exactly (actually SUg ~ 1 eV). In intentionally asymmetric quantum wires and 
dots with strong confinement, this cancellation will be mitigated, and more pronounced 
spin splittings will be attained. 



4-2. Wurtzite semiconductors 

For Brillouin zone-centre conduction band states in structures based on wurtzite 
semiconductors, the band edge function Uso transforms as belonging in the Fi 
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representation of the space group Cgt,. It is formed of atomic s and p orbitals, so 
that and Hf°'^ are finite |12]. 

The basis functions m„o can still be chosen real, producing the identities ( IT6l) . Then 
we have, for the 3rd perturbative order from (fT2!) . 



(25) 



We have used the identity 

kaki3 f - f kakfs = -ika (V/j/) - i (Vaf) kp. (26) 

For the fourth perturbative order, similar to the zinc-blende case, the second term 
on the right-hand side of equation (ITS!) gives no contribution. We have 

Hif =B^^p [VJV kp - VpW ko\ + 52„^ [k^ {VpW) + {VJV) kp] 
+ [VJ kp - Vpf k^] + C2ap [kc (V^f) + (VJ) kp] . 



(27) 



where 



-Dla/3 — > : 2 2 ' '■ 2 2 ' \'^°) 

B2ap = — —277—2 ' y^^) 



^lap — / ^ : 2 ' '■ 2 

zniQ Usi Usv cJsi" ^rn^ ^si ^sV ^si" 



(30) 



^ _ Sr^ ^{Pa) slip 15)11' SUi'i'f H^,^^ , h^H^i \pa)ll' {Pf})l'l" SUl"s .„^s 
<-^2a/3 — / ^ : 2 ' '■ 2 ' y'^^) 

j^^^ im^Q Usi uJsi' ojsi" Usi tOsV ojsi" 

Finally, letting the wurtzite c-axis be along the z-direction, we have Hgo = 

£>(3) , fV(4). 
J^ss \ ss • 

Hso =Riwz ([V/ X k]^. or, + [V/ X k]^ ay^ + Rl^^z [V/ x k]^ a, 

+ ai [(V,/) k, + ky (V./)] - «i [(V./) k, + k^ (V J)] a, 

+ «2 [(V,/) ky + A;, (V,/)] (T, ^^^^ 

+ i?2W^z ([VW^ X k]^ a, + [Viy X k]^ a,) + i?;^^ [Viy X k]^ a, 

+ /3i [{VyW) k, + A;, (V.l^)] a., - /3i [(V.iy) A;, + k^ {V ,W)] ay 

+ /32 [(V.iy) + (V,iy)](T,. 

The material parameters entering here are given in appendix B. The structure of this 
complicated expression resembles the net SIA spin-orbit Hamiltonian for zinc-blende 
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systems (1191) . The differences are plain. Due to the anisotropy of the wurtzite, 
Riwz 7^ R'lwz R2WZ 7^ R^wz Rashba terms. Another difference consists 

in the presence of new contributions due to the finite matrix elements of the spin-orbit 
interaction H^^'^"^ and H^fi'^"^ between the conduction and remote bands. New terms are 
proportional to ai, a2, Pi and /32, with ai 7^ 0:2, /3i 7^ 1^2 due to the anisotropy of the 
wurtzite. These terms are reduced to the conventional form of the Rashba spin-orbit 
interaction for ID external electric field and variation in the chemical composition, if 
we put VxW = VyW = Vxf = Vyf = 0, VzW ^ and V2/ 7^ 0, which is easily seen 
from expression (1321) . 

It is interesting to learn which bands contribute to make ai, 02, /3i and finite 
(see appendix B). It can be deduced using the tables of direct products of irreducible 
representations of the space group Ceu, see [47j: Fg x Fg = Fi -|- r2 + Fg, F2 x Fg = Fg, 
and Fi x V j = Tj for any j. We should also take into account that polar vectors (e.g. 
k) transform as Fi + Fg, while axial vectors (e.g. [Vt/i x k] ) transform as F2 + Fg 
[9]. Then we immediately conclude that bands with symmetries Fi and Fg define the 
strength of the parameters ai and Pi. Hence, they are finite even in the truncated 
bands Kane-like model with the nearest valence bands Fg and F^ The term involves 
new matrix elements, not expressed via known band parameters [M]: (Fg | hky \ F^) and 
(s I [VUi X k]^. I F^), see also [12]. 

For parameters 02 and /92, in a similar way, there should be finite contributions from 
the states with symmetries F2 and Fg. This turns ^2 and P2 to zero in the truncated 
bands model. The bands F2 do not appear in the pseudopotential calculations [18| H9]. 
which probably means that they are very remote. 

Let us consider only the nearest valence bands Fg and F^ to estimate Riwz, R'iwzj 
R2WZ and R'2wz GaN/AlN. We can use the expressions fl23l) and flMl) . As previously, 
the material parameters are taken from [31]. We obtain R2WZ = R'2wz ~ —0.01 A^. 
Accidentally, due to different signs and comparable strength of the terms entering 
the parenthesis of expression (!23|) . we have vanishingly small Riwz = R^iwz ^ 
0.7 meVA^. The external electric field-induced spin orbit interaction dominates because 
the conduction band offset SUg ~ 2 eV, so that \R2wz5Us\ = 20 meVA^ > Riwz- While 
the value of the effective Rashba parameter for GaN/AlN is only one- hundredth of that 
for the GaAS/AlAs system, very strong electric fields acting on electrons in GaN/AlN 
can induce large spin splitting of electron states, comparable to that in narrow-bandgap 
materials |28j . 

5. Conclusions 

In the k ■ p method, we derived SIA spin-orbit interaction terms for conduction 
band states near the Brillouin zone centre in zinc-blende and wurtzite semiconductor 
heterostructures taking into account all remote bands. The results are applicable 
to quantum wells, wires or dots. Electric field-induced terms and heterointerface 
contributions were considered, both generally having comparable strength. They can 
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be written in a unified manner only for 2D systems. Tlie resulting expression for the 
spin-orbit Hamiltonian f|T9|) in zinc-blende materials takes the conventional form of 
the relativistic spin-orbit interaction [1]. The Kane model is adequate to establish this 
form, with other remote bands only changing the values of the material parameters. For 
wurtzite materials, the net SIA spin-orbit Hamiltonian (|32l) has a complicated form due 
to the anisotropy of the wurtzite and new contributions, which appear owing to finite 
matrix elements of the spin-orbit interaction H^^'^°^ and -ffff**"^ between the conduction 
and remote bands. Knowledge of the parameters of the Kane model alone is insufficient 
to write the spin-orbit Hamiltonian for wurtzite. The effect of remote bands is yet to 
be evaluated. 

We analysed two popular semiconductor pairs, GaAs/AlAs and GaN/AlN, with 
the goal of establishing the mechanisms actually governing SIA spin-orbit interaction 
in heterostructures composed of these materials. Both pairs have 'accidental' sets of 
parameters strongly differing from the 'general' picture. For 2D systems GaAs/AlAs, 
the interface and external electric field-induced contributions are comparable and have 
different signs partially canceling each other and significantly reducing the net spin 
splitting. In asymmetric quantum wires and dots with strong enough confinement, 
this cancellation will be mitigated, and relatively more pronounced spin splittings can 
be attained. For the systems based on GaN/AlN, the nominally interface-induced 
contribution is very small as compared to the external electric field-induced one. The 
evaluations were based on parameters available for the Kane model alone [31], without 
remote bands. 
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Appendix A. Lowdin perturbation scheme up to the fourth order 

We sketch the Lowdin perturbation scheme [33], which is an efficient tool for the k ■ p 
diagonalization treatment [TH |35l SSI |16]. We follow [9] and then derive all necessary 
elements for the fourth-order correction. The k • p system ([8]) can be presented as: 



where 6nn' is the Kronecker delta, and H' = Hi -|- H2 is the perturbation. The all-band 
Hamiltonian is being decomposed into the block of m-indexed states, whose mutual 
interaction is taken into account exactly, and a block of 'remote' /-indexed bands, treated 
as a perturbation. Here our m-class block consists of only a single conduction band 




(A.l) 



Here Hq is the Hamiltonian of the zero-order approach: 




(A.2) 
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m = s, but we are preserving the general notations with indices m, m', m", etc all 
belonging to the m class, following [9] . The perturbation Hi does not contain elements 
of interaction of m-indexed bands with other bands. In other words, it is already block- 
diagonal. It is these elements of m-l band interaction that are contained in H2. 
The canonical transformation of the envelope functions 

A = e-^A (A.3) 

with an anti-Hermitian (S"*" = — S) matrix S results in a set of equations 

H A = eA (A.4) 

with 

H =e-^He^. (A.5) 

The properly chosen S must provide a decomposition of the whole system into a set of 
equations for m-bands and an abandoned set of equations for other bands. Expanding 
exp(S) in a series 

eS = l + S + ls2 + is3 + ..., (A.6) 



we obtain 



Here 



00 

H = 5^-{H,S}("). (A.7) 



n 

n=0 



{H,S}(°)=H, {H,S}(^) = {H,S}, (A.8) 

{H,S}(^) = {{H,S},S},... (A.9) 
Let us consider 

8 = 81 + 82 + 83, (A.IO) 

where 8„ is the matrix of the nth order in H'. The elements 8„ are defined using the 
recursive equations 

{Ho,Si} + H2 = 0, (A.ll) 
{Ho, 82} + {Hi, 81} = 0, (A.12) 

{Ho, S3} + {Hi, 82} + \ {{H2, Si} , Si} = 0. (A13) 

Then the transformed Hamiltonian H takes the form 

H =Ho + Hi + ^ {H2, Si} + ]- {H2, S2} 

2 2 ^^^^^^ 

+ -{H2,S3}--{{{H2,Si},Si},8i} 

to the fourth order in H' inclusively, having no elements of m-l band interaction. It is 
interesting to note that an apparently similar perturbation method given by Luttinger 
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and Kohn [30] actually differs from the Lowdin's in an attempt to make a redundant 
diagonalization inside the remote bands block. 

The matrices Si and S2, which are derived from equations ( lA.llI) and ( 1A.12I) . are 
known [9]: 

Simi = (A.15) 
where Unn' = Cno - Cn'o, and 

m' (' 

Now using equations ( 1A.13I) . ( 1A.15I) and (1A.16I) . along with the anti-Hermiticity of 
S, we obtain 

^■imi — y — 5 \ H 



m'V 



+ y^EL!A2iEi( — I — + — I — ] .^17) 

TTl TTl TJt TT/ TT I TT I 

^rnm'-'^rn'm"-'^m"l \ ^ ^ ml' -"^ I' I" -"^ I" I 



_ -'-'-mm'-'-'-m'm"-'-'m"< _ 

m ,m" I' ,1" 

Appendix B. Material parameters entering expression ( I32|) 



rilWZ — > A- A 2 



1,1' - '0 



-At 



^ ^ {Vy)sl ([V^l X k]J^^, bUyy, {p,)^,. 



fi\py)siSUu' {p.)i'i" ([V^ixk] 



xn"s 



p/ _^ ^ {PX)SI {[^^U X k]j„, {Py)vs 

.2,,,,,,, „ 



_^ K^{Px)si5Uii' {[VUi X k]^)^,^„ {py)i"s 

2im^ c^ujsi ujsi' Usi" ^^''^^ 



+ E 



([Vf/i X k] J^, {p^)u, 5Uin" {py)i"s 
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^h\py)sl{p.)u' ([V5C/Xk]J,^ 



+ 



+ 



h'{[V5Ux\,]^)^^ {py)u'{p.)vs 

Aiml c2 Usi ujsi' 
{[WUi X k]^)^^ {py)w {pz)i'i" SUi»s 

Aimf^ to si CO si' ^si" 
h'^SUsi {py)ii' {pz)i'i" {[^Ui X k]^)^„^ 
Aiml Usi Usi' ujsi" 

h'^{Py)sl {Pz)ll' SUi'i" ([V^l X k]^)^„^ 

4im^ Usi ujsi' ^si" 
hf" ([Vt/i X k] J^, 6Uiv {py)vi" {Pz)l"s 



Oi2 



= E 



'^^ ^ti 

nmvu.x ki)^^ {p.)w{py)i 



+ 



ft = E 



h\Py)si{p.)w ([V^ixk]j,^ 
h'{[VU,xk\^)^^{py)u>{p;)i 



(B.3) 



4im^ c2 LJsi" 

h\p.)si{Vy)w {\V5Ux^;)^,^ 
2iml Usi LJsi' 

^ ([V^l X k] J^^ (p^);;/ SUi»s (3 4) 

2im^ c2 Usi' to si" 

^ h'^{Px)sl {Py)lV SUlH" {[VUl X k] J^„^ 
2imo tU^i tUs;/ (J^;// 

P ^ fiHPy)si ([V^i X k]J^^, f 1 1 

^2WZ — / ^ -r. I —2 \ 2~ 

4im^ c2 c<;2 
4im^ c2 ujsi f^si' 

^\Px)sl {Py)lV ([V^l X k] J^,^ ^g_g^ 



(B.7) 



Spin-orbit interaction in three- dimensionally bounded semiconductor nano structures 17 



E 

1,1' 



(B.8) 
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